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I. INTRODUCTION 



At an early stage in the development of quantum mechanics, Schrodinger at- 
tempted to find nonspreading wave-packet solutions for a quantum-mechanical parti- 
cle evolving along a classical trajectory. He succeeded for the case of a particle subject 
to a harmonic-oscillator potential, obtaining a solution that follows the classical os- 
cillatory motion without changing shape M. This solution is now called a coherent 
state 0. Schrodinger also attempted to find analogous solutions for the Coulomb 
potential, without success B. Many authors since have discussed this issue, and 
it is now known that there are no exact coherent states for the Coulomb problem 
§ |, |, 0, §,§. 

Although the original Schrodinger problem for the Coulomb potential has no solu- 
tion, one can nonetheless obtain minimum-uncertainty wave packets exhibiting many 
features of the corresponding classical motion. For example, it has recently been 
shown that a type of squeezed state called an elliptical squeezed state (ESS) is a 
minimum-uncertainty wave packet for the planar Coulomb problem fllPf . An ESS is 
localized in both the radial and angular coordinates, and it travels along a keplerian 
ellipse in two dimensions. The radial part of the initial solution is a radial squeezed 
state (RSS), which minimizes the uncertainty relation for a set of radial variables 



that express the Coulomb problem in a form similar to that of an oscillator HIT], 12 



The angular part of the initial solution is a circular squeezed state (CSS), which 
minimizes the uncertainty relation between the angular momentum L and a suitable 



angular-coordinate operator [10], [13 



The primary goal of this paper is to extend the ESS construction to three dimen- 
sions and thereby to provide a class of wave packets coming as close as possible to a 
solution of the original Schrodinger problem for the Coulomb potential. We call these 
wave packets keplerian squeezed states (KSS). We show that the KSS are minimum- 
uncertainty wave packets that are localized in all three dimensions and that travel 
along a keplerian ellipse. Three-dimensional wave packets of this kind are of partic- 
ular interest at present because experiments using short-pulsed lasers are attempting 
to produce Rydberg wave packets that move along elliptical orbits. 
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To date, experiments that have detected the classical oscillation of a wave packet in 
a Coulomb potential have been performed using purely radial wave packets consisting 



of a superposition of n states with I fixed, typically, to a p state WQ ITR [161 II7L [T3L [19 



These states follow the initial classical motion, but also exhibit distinctive quantum- 
mechanical features. For example, the uncertainty product oscillates periodically as a 
function of time. This is a characteristic of squeezed states [fZlJ, and indeed the RSS 
provide a good analytical description of wave packets of this type. The wave packets 
also exhibit quantum interference effects, as they undergo a series of collapses and 
full/fractional revivals and superrevivals fl2~T|, |22], ^3], ^4], [26[ . 

To generate a three-dimensional wave packet localized in radial and both angular 
coordinates, a superposition of n, I, and m levels must be created. This requires the 
presence of additional fields during the excitation process. One proposal for achieving 
this involves using a short electric pulse to convert an angular state into a localized 
Rydberg wave packet moving on a circular orbit |2"?fl . An additional weak electric 
field could then distort the orbit into an ellipse of arbitrary eccentricity. We expect 
the motion of these wave packets to be well described by the KSS we present here. 

The construction of the KSS requires the identification of angular operators ap- 
propriate to a problem with spherical symmetry. These can then be used to obtain a 
class of squeezed states, called spherical squeezed states (SSS), that minimize uncer- 
tainty products for the angular variables. This construction and the basic properties 
of the SSS are presented in Sec. II. 

Section III gives the construction of the KSS. At the initial time, the KSS may 
be written as products of RSS with SSS. The resulting wave packet is localized in all 
three dimensions and consists of a superposition over the three quantum numbers n, 
I, and m. In this section, we also calculate expectation values for physical quantities 
characterizing the initial wave packet, and we discuss its time evolution. 

Section IV outlines how the KSS can be generalized to include the effects of 
quantum defects. This is necessary because experiments performed on Rydberg wave 
packets generally use alkali-metal atoms rather than hydrogen. 

Section V contains a summary of our results. Some mathematical formulae used 
in deriving expectation values in Sec. II are given in the Appendix. Throughout this 
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paper, we use atomic units with U — m — e — 1. 



II. SPHERICAL SQUEEZED STATES 

In this section, we obtain the spherical squeezed states. These are squeezed states 
for angular coordinates on a unit sphere. Section IIA discusses the choice of angular 
variables used in the construction. Section IIB gives the derivation of the SSS, and 
Sec. IIC examines their properties. 

A. Quantum Angular Variables on the Sphere 

Classically, the standard angular coordinates on the unit sphere are the spherical 
polar coordinates 9 and 0, which can be taken as defined in the intervals < 9 < tt 
and —it < < 7r, respectively. At the quantum level, however, 9 are are not well 
suited as coordinate operators. Some of the difficulties stem from problems already 
appearing in the planar case. For example, if is taken as periodic, the discontinuity 
at = it causes complications when is acted upon by derivative operators such 
as the z component L3 of the angular momentum L. The alternative assumption 
of a continuous angle in the interval —00 < < 00 also introduces complications 
because this variable is not periodic and hence is not observable in the Hilbert space 
for which L3 = —id^ is hermitian. In the present instance, the use of in conjunction 
with 9 on the unit sphere raises further complexities because is not well defined at 
the poles 9 = 0, it. 

Difficulties with angular coordinates on the circle have been widely discussed in 



the literature po |. In addressing the Schrodinger problem for the planar Coulomb 



potential [|Hj and in a related work on minimum-uncertainty angular packets [[13 



we followed ref. [£9] in circumventing these and other difficulties by replacing with 
two quantum coordinate operators, sin0 and cos0, which are both continuous and 
periodic. Use of these variables permitted a completely satisfactory resolution to the 
essential issues. 

In the present work, we follow a similar strategy. We introduce three angular- 
coordinate operators a,j, j = 1,2,3, via their matrix elements obtained by insertion 
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of the three continuous and periodic coordinate functions fj = (sin 9 cos </>, sin 9 sin 0, 
cos 9), respectively, in the Hilbert-space inner product: 



where dfl is an element of solid angle on the sphere. They obey the relation J2j a "j — 1- 
Classically, they are a natural choice corresponding to the identification (x, y, z) — > 
(sin 9 cos <j), sin 9 sin <ft, cos 9), which uniquely covers the unit sphere. Moreover, this 
choice reduces in the limit 9 — > 7r/2 to the planar case. 

The quantum operators dj and the components of the angular-momentum 
operator L obey the commutation relations 



Each of these commutators has an associated uncertainty relation that must be con- 
sidered in the search for minimum-uncertainty solutions. 

The number of nontrivial relations can be substantially reduced by stipulating 
certain initial conditions on the desired wave packet. Without loss of generality, 
we take the initial minimum-uncertainty configuration to be a wave packet localized 
about the positive x axis. Since we seek a packet moving along an orbit in the x-y 
plane, we also assume the initial shape is reflection-symmetric about this plane. It is 
also physically reasonable to suppose that the initial packet is reflection symmetric 
about the x-z plane. These requirements impose the coordinate conditions 





[dj, Lk) — iejkiai 



(2) 



(a 2 ) = (a 3 ) = 



(oi) > 



(3) 



and the angular-momentum conditions 



(Li) = (L 2 ) 







(4) 




and 




(6) 
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B. Construction of SSS 



The CSS are defined as states minimizing the relation (pj) in the x-y plane QlO 



Requiring that the SSS reduce to the CSS in the planar limit ^ — ^ 7r/2 means that 
(H) must be minimized in three dimensions too. Our construction of the SSS in 
this section therefore begins with this minimization, after which we return to the 
remaining relation (|6]). 

A wave function X {9, 0) minimizing Eq. (|5|) obeys the differential equation 

a 2X (9A) = -i](L3-(L 3 ))x{0,<f>) , (7) 

where the squeezing 5 in the angular coordinates is a real constant given by 

1 _ A« 2 _ 2(A« 2 ) 2 

5 AL 3 ~ !(«!>! " • 1 ] 

As expected, setting 9 — > it/ 2 reduces Eqs. (|7|) and (g) to the defining equations for 
the CSS. Fixing (L 3 ) = /3, a real parameter, Eq. (0) becomes 



+ 5sin#sin0-i/3)x(#,0) = . (9) 

The solution to this equation is 

x (9,<f>) = N exp (5 sin 6 cos (f) + i{3(f) + f (6)) , (10) 

where f(9) is an undetermined complex function of sin# and cos^. For x{9, <p) to be 
single valued, f3 must be an integer. This follows from the general result |13] that 
minimum-uncertainty angular wave packets must have integer expectation values (L3) 
of the angular momentum. 

We next address the issue of the role of the second uncertainty relation (0). An 
immediate question is whether the function f{9) may be chosen so that x{9,4>) also 
minimizes (||). This would require that, in addition to @, x{9,4>) must satisfy the 
differential equation 

cos 6* x(9, 4>) — *T7 ( — i cosipde + % cot 9 sin <p dj,) x(9, 4>) , (11) 


where 5' is a second squeezing parameter defined by 1/5' = Aa 3 /AL2- However, 
substituting fllOD into ( pTT|) produces an equality only if 5 = 5', (3 — 0, and f{9) 
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is a constant. Physically, this means the squeezing ratios are equal, Aa 2 /AL 3 = 
Aa 3 /AL 2 , and that the mean angular momentum vanishes, (Li) = (L 2 ) = (L 3 ) = 0. 
However, this is unacceptable because we seek solutions with a nonzero value for (L 3 ). 
Indeed, the two uncertainty relations (|) and (|6|) can be satisfied simultaneously only 
for states symmetric under rotations about the x axis. The rotational symmetry must 
therefore be broken, and the relation (|5]) cannot be minimized. 

One simple way to break the rotational symmetry would be to keep (L 3 ) ^ 
in Eq. fllPl) but to set f{9) = 0. The function x(@,4>) then would take a simplified 
form that reduces directly to a CSS as 9 — > tt/2. This choice has the added elegant 
feature that expectation values would take the same form as those for the CSS given 



in ref. ||10|| , except with modified Bessel functions l v replaced by modified spherical 
Bessel functions i v with different arguments. However, wave packets of this type are 
also unacceptable because certain physical expectation values diverge. For example, 
(L 2 — Lg ) diverges. 

Obtaining an acceptable minimum-uncertainty wave packet evidently requires set- 
ting (3 7^ and deriving a suitable nonconstant f(9). Further physical input is needed 
to accomplish this. Note that, whatever its form, the function f(9) must satisfy the 
limit f{9) — > as 9 — > tt/2, so that x{®-> 0) reduces correctly to a CSS. Furthermore, 
the reflection symmetry about the x-y plane suggests that / = /(sin#), independent 
of cos 9. 

The function f(9) controls the spread transverse to the orbit of the initial wave 
packet. This spread should depend on the value of (3 = (L 3 ), since physically as (L 3 ) 
increases the wave packet is expected to become more confined to a region around 
the x-y plane. For simplicity, we take the spread as independent of the spread around 
the orbit, which is controlled by the parameter 5. It turns out that this assumption 
suffices for our purposes. If it is relaxed, the resulting wave packets have a significantly 
more complicated analytical structure. In what follows, we therefore permit f(9) to 
depend on /3, but we require it to be independent of 5. 

A wave packet following a keplerian orbit with maximal confinement to the x-y 
plane should have minimal expectation value \(L 2 — L 3 )\. We can use this physical 
condition to obtain f(9), as follows. Since f{9) is independent of 5, the limit <5 — > 
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can be imposed in Eq. (|T0|) without loss of generality. In this limit, the uncertainty 
AL 3 vanishes because x{@-> 0) becomes an eigenstate of L 3 with eigenvalue m = f3. 
Expanding 4>) m terms of spherical harmonics Y[ m (9, (j>) with m = (3 gives 



= E c iVM) • (12) 

5^0 Z>/3 



shows that 

5^0 



Then, calculating the expectation value | (L 2 — L\) | in the state 
it is minimized when the coefficients q in Eq. ( |12|) are proportional to Sip. The sum 
reduces to a single term proportional to Ypp and hence proportional to skr 9 exp(i/30). 
This implies that exp(/(#)) = sin' 3 9. As expected, f(9) is a function of sin# and 
satisfies the condition f(n/2) = 0. Moreover, this particular functional form leads to 
an SSS without singularities at the poles and hence with finite physical expectation 
values. 

Combining this result with Eq. (|Tt]) produces a set of minimum-uncertainty states 
satisfying the requirements for SSS. We therefore take the SSS to be given by 

X{0,4>) = iVsin /3 flexp(5sinflcos0 + i(3<p) . (13) 

C. Features of SSS 

We next calculate the normalization constant and some expectation values for the 
SSS. These are most conveniently specified in terms of a function Aj(5), defined as 

A?(6) = 2tt I" dO (sin 9) w+]+1 Ij(25 sin 9) , (14) 

w 

where j is an integer and Ij{z) is a modified Bessel function of the first kind. The 
Appendix presents analytical expressions for the functions A? (5) and discusses some 
useful properties. 

The normalization constant iV is given by 

N = ^r= ■ ^ 

The nonvanishing expectation values of the angular-coordinate operators are 
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(*l) = ^f7rA^ +1 (5)+A$(5)) , (17) 
2A (d) 



2x 1 A{{5) 1 



{al) = T5 4^ = Ts {s ' mdcos(p) ' (18) 

(al) = -^-jA%(6)-A% +1 (6)) , (19) 
A o{d) 

while some useful expectations involving the angular momentum are 

< L 3> = !^ + /^<^os0>+/3 2 , (20) 
2A {b) 2 

(L 2 ) = (3((3 + 1) - 5 2 (l- ((sin^ cos 0) 2 >) 

+ + . (21) 

Using the identity ( |55| ) for the (<5) in the Appendix, one can show that J2j a<2 ' 1 



2\ 



as required. Note that, in the limit 5 — ^ 0, (L 2 ) = (3 and (L ) — + 1), as is 
appropriate for an eigenstate with I = m = (3. 

We next turn to a discussion of some properties of the SSS. First, expanding x(#, (f>) 
in spherical harmonics for the case with 5^0 shows that it is a superposition of both 
/ and m states. The condition exp(/(6 1 )) = sim ensures that the squared modulus 
\ci m \ 2 of the weighting coefficients c\ m in the expansion is maximal for / = m = (3 
and falls rapidly away from these values. The rapid decrease indicates that (L 2 — L 2 ) 
remains small for 5 ^ 0, as expected. The reader interested in the form of the 



expansion can find it in Eq. ( |3T| ) of Sec. III. 

We have seen that the parameter (3 is the expectation of the z component of the 
angular momentum. It turns out that the parameter 5 has an interpretation similar 
to that of its counterpart for the CSS. The uncertainty AL3 depends on 5. This 
may be confirmed in the neighborhood of x — 1, y = z = on the unit sphere 
by expanding the wave packet around <p = at 6 = ir/2. We find |x(7r/2,0)| 2 oc 
exp[25(l — \4> 2 + - ■ •)], which to leading order is a gaussian in (ft with standard deviation 
a = 1/V26. Therefore, for larger values of 5 the initial wave packet becomes narrower 
in <b. 
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Classically, a particle moving along a trajectory confined to the x-y plane has 
angular motion depending only on 0. The trajectory is initialized by giving values 
of sin0, the sign of cos0, and L 3 . At the quantum level, the initial angular position 
of an SSS is specified by the expectation values in Eqs. (|3|). As we have seen, the 
parameter j3 gives the initial angular momentum (L3) and controls the transverse 
angular spread, while 6 controls the angular spread along the orbit. Compared to 
the corresponding classical problem, a general quantum packet would need two ad- 
ditional parameters determining the spread on the surface of the sphere. However, 
our simplifying assumption for f(6) means that only one extra parameter is needed 
to specify the quantum solution for the SSS. A related point is that the SSS and CSS 
depend on the same number of parameters. This is physically reasonable since both 
states follow the same keplerian orbit in the two-dimensional x-y plane. 

III. KEPLERIAN SQUEEZED STATES 

This section discusses the keplerian squeezed states. Their definition is presented 
in Sec. IIIA, where their parameters are determined in terms of specified physical 
quantities. The time evolution of a KSS is examined in Sec. IIIB, and an example of 
a KSS wave packet moving along an elliptical orbit is provided. 

A. Construction and Specification of KSS 

The SSS solution flllf ) is a function only of the angular coordinates and has neither 
time dependence nor dependence on the quantum number n. A solution consisting of 
a product of an SSS with a radial energy eigenstate of given n might be considered, 
but this produces a stationary state. However, we can create a three-dimensional 
state moving on a keplerian orbit by combining an SSS with an RSS. 

The RSS are constructed and analyzed in refs. [[II], They are wave packets 
localized in the radial coordinates that initially undergo oscillatory motion between 
the inner and outer apsidal points of the corresponding keplerian ellipse. The pe- 
riod T c \ of the motion is that of a classical particle moving in a Coulomb potential. 
The construction involves converting the classical effective radial hamiltonian for the 
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Coulomb potential in terms of conventional radial variables to an oscillator descrip- 
tion in terms of new variables. The resulting classical problem is quantized, and wave 
functions are obtained that minimize the ensuing quantum uncertainty relation. The 
RSS are given by 

iP(r) = N'r a exp[- lo r - i lir ] , N' = [(2 7o ) 2a+3 /r(2a + 3)] 1/2 . (22) 

For the initial wave packet, the parameters a, 70, and 71 determine the radial position, 
the uncertainty in the radial variables, and the radial momentum. In refs. |JTl~l , [12| , 



the RSS are shown to describe a Rydberg wave packet that has been excited by a 
single short laser pulse. The angular part of the full three-dimensional wave function 
in this case is fixed to a p state for excitation from the ground state. 

A minimum-uncertainty wave packet moving on a keplerian orbit can be obtained 
by combining an RSS with an SSS. The resulting wave packet is a KSS. The sepa- 
rability of the full hamiltonian and the independence of I of the RSS uncertainty in 
the new radial coordinate together make it possible to minimize simultaneously the 
uncertainty relations for the new radial variables and those for the angular variables 
in (|5|). We therefore can take as an initial wave function the product of ip(r) in ( p2[ ) 
and <p) in flTB|). The result is a normalized five-parameter family of KSS, 



= NN'r a sin p 6 exp [5 sin cos - ( 7o + z 7l )r + , (23) 

where N is given by Eq. flT5|) and N' is given by Eq. (^2|). The KSS are minimum- 
uncertainty wave packets localized in all three dimensions. The choice of the initial 
angular-coordinate location is implicit in the SSS construction and is specified in Eq. 

(!)• 

Expectation values of operators for the KSS can be calculated analytically. Since 
the radial and angular wave functions separate, the angular operators have the expec- 
tation values given in Eqs. ( |T6| ) - fl2~T|). The expectation values for the radial operators 
are 

W = ^ , <i> = ^T , (24) 
270 r a + 1 
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(a + 2)(2a + 3) 
2V 




2 7 o 2 



(25) 



(a + l)(2a + l) 



(Pr) 



71 , (Pr 2 ) 



2a + 1 




(26) 



The uncertainty for the radial coordinates is 



ArAp r = \ 



2a + 3 
2a + 1 



(27) 



The RSS are not minimum-uncertainty states in the variables r and p r , which is as 
expected since the construction minimizes the uncertainty in the new radial variables 
instead. For large values of a, however, ArAp r —> 1/2. 

The expectation value for the energy (H) is obtained using the full Coulomb 
hamiltonian and depends on all five of the KSS parameters. It is 



We next address the issue of initializing a KSS. With the wave packet located at the 
outer apsidal point of an elliptical orbit, the uncertainty product for a radial Rydberg 
wave packet is a minimum. Initially imposing the constraints (||) ensures that the 
ellipse has semimajor axis aligned along the x axis. The five KSS parameters a, /3, 70, 
71, and 5 can then be fixed by specifying the expectations of the radial coordinates 
(r) and (p r ), the expectation of the angular momentum (L3), the expectation of the 
energy (H), and the spread AL3 in the angular momentum. 

The natural choice for the initial expectation value of p r is zero. Similarly, the 
natural choice for the initial expectation value of r is the outer apsidal point of 



the orbit, r out = n 2 [1 + y (1 — /(/ + l)/n 2 j. The expectation (L 3 ) gives the initial 
angular momentum of the wave packet and its spread transverse to the orbit. The two 
remaining conditions fix the initial spread of the wave packet in the radial and orbital 
directions. For the first, a natural choice is to set (H) equal to the mean energy of a 
Rydberg wave packet consisting of a superposition of n states centered on the value 
n. A packet of this type is produced by excitation with a short laser pulse tuned to 
the mean energy E n = —l/2n 2 . For the second, it suffices to specify the width AL 3 
of the superposition. Note, however, that a single laser pulse cannot excite a wave 



(H)= 1 -( P 2 r ) + l(- 2 )(L 2 )-(-) 
2 2 r l r 
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packet localized in the angular coordinates. Creation of a superposition of angular 
eigenstates requires additional fields to mix angular-momentum eigenstates. 
The full set of conditions sufficient to fix the five KSS parameters are then 

(r) = r out , (p r ) = , (L 3 ) = (3 , 

(H) = E n , y/(Ll) - (L 3 ) 2 = AL 3 . (29) 

These determine the values of a, (3, 70, 71, and 5 in terms of physical parameters in 
the excitation process. 

B. Evolution of KSS 

We next investigate the time evolution of a KSS wave function matched to a 
Rydberg wave packet at the outer apsidal point of an elliptical orbit, as described 
above. By construction, the packet is expected to travel along a classical ellipse in 
the x-y plane. 

Since (p r ) = and (L 3 ) = (3, the initial motion of the wave packet for (3 > is 
in the direction of increasing 0. The geometry of the ensuing orbit depends on the 
values of n, (3, and AL3. The parameter (3 gives the central values of / and m. For 
(3 ~ n — 1, the KSS orbit becomes circular, with the wave packet propagating at fixed 
mean radial distance from the origin. As (3 decreases, the orbit becomes elliptical, 
with the inner apsidal point moving closer to the nucleus. Also, for /3 — > 1 the radial 
wave function becomes oscillatory as the electron passes close to the nucleus. 

To study the time evolution of a KSS, we expand $?(r,9,<fi) in (p3|) in terms of 
energy and angular-momentum eigenstates, 

V(r,9,<j>,t)= CnimRni(r)Y lrn (e,(f))e- iEnt . (30) 

n,l,m 

The expansion coefficients c n i m = (&(r, 9, <fi, O)|i2 n i(r)Yj m (0, </>)) can be calculated us- 
ing Eq. (pf) as the initial wave function ty(r,9,(f),0). The radial and angular parts 
separate, and we may write c n \ m = c£j , where = (tjj(r)\R nl (r)} for ip(r) in 
(H) and c£ s) = (x(M)l^m(M)) for x (9, (/)) in (0). The coefficients c^ ss) for the 
radial part of the expansion are given in Eq. (58) of ref. [0. The coefficients for the 
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angular expansion can be calculated using The result for m > is 



21 + 1(1- m)\ 
\ An (l + m)\ 



m) min(m,/3) 





N^5 , { l-m) mod2 (-1) 



E E 



2 k k\ 



k=0 p=0 



X 




p ^min(m,/3)^ ^ 



z i 



(31) 



Here, A/" is the SSS normalization constant in (|T^), min(m,/3) gives the minimum 
of the two values m and (3, and i v (z) is a modified spherical Bessel function. The 
appearance of the Kronecker delta <5o,(«-m) mod2 implies that if I — m is an odd integer 
then cf^ = 0. Therefore, the SSS is composed only of eigenstates for which I — m 
is even. This follows because the associated Legendre functions with odd I — m are 
odd functions of 9 in the interval < 9 < 7r centered on n/2, while <fi) is an even 
function of 9 in this interval. 

As an explicit example of the evolution of a KSS, consider matching \l/(r, 9, 0, 0) 
to a Rydberg wave packet at the outer apsidal point with n = 45, (L3) = 30, and 
AL3 = 2.5. Using these values and Eqs. (p9|) , we obtain the KSS-parameter values 



a ~ 62.846, [3 = 30, 70 ~ 0.01834, 7l = 0, and 5 ~ 12.826. This gives (r) = 



number of terms with n centered on n and I and m centered on (3. In the present 
example, we allow a spread of 10 units in n and / and four units in m. This gives 484 
coefficients, half of which vanish. We keep the remaining 242 terms in the series and 
plot two-dimensional sections through the KSS as a function of t. 

Figure 1 presents a slice through the initial KSS in the x-z plane. This slice is 
transverse to the classical orbit, which lies in the x-y plane. The figure shows that 
the initial wave packet is localized around a point on the x axis at r w r ont and in a 
narrow range of 9 near 9 = tt/2, corresponding to z ~ 0. 

Figure 2a presents a slice through the initial wave packet in the plane of the 
classical orbit. The initial wave packet is located on the positive x axis at the outer 



r out ~ 3508.6 a.u. and (L 2 ) ~ 938.1. Defining a mean value / of / from the relation 
(L 2 ) = 1(1 + 1), we find f ~ 30.1 w (3. 

The series in ([30]) may be well approximated by truncating the sum to a finite 
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apsidal point. Figures 1 and 2a taken together show that the initial KSS is localized 
in all three dimensions. 

The classical keplerian orbit for a particle in a Coulomb potential is T c \ = 2irn 3 . 
With n = 45, we obtain T c \ ~ 13.4 psec. Figure 2b shows the KSS in the x-y plane 
at t — T c i/3. It has moved in the direction of positive <ft an d is spreading along the 
elliptical orbit. In accordance with Kepler's laws, the wave packet moves more slowly 
near the outer apsidal point than near the inner one. As a result, it has traveled less 
than 1/3 of the orbital circumference at t — T c \/2>. 

Figure 2c shows the KSS in the x-y plane at t = T c \/2. It has spread along the 
elliptical orbit and is moving more rapidly. The radial distance to the inner apsidal 
point is r in ~ 536 a.u. This is sufficiently far from the nucleus for the wave packet to 
remain localized in r, and hence no radial oscillations are apparent. 

Figures 2d and 2e show the slice through the KSS in the x-y plane at the times 
t = 2T c \/3 and t = T c \, respectively. The motion slows again as the wave packet 
approaches the outer apsidal point and becomes more localized. At t = T c \, the wave 
packet resembles the initial wave packet. However, the motion is not exactly periodic. 
As time increases, the wave packet collapses and for t > T c i a cycle of revivals and 
superrevivals commences. 

Figure 3 shows the wave packet in the x-y plane at t = T c \/2, but viewed from a 
point on the positive x axis looking towards the nucleus. The elliptical shape of the 
orbit is evident. 



IV. INCORPORATION OF QUANTUM DEFECTS 

Experiments studying the behavior of Rydberg wave packets are usually performed 
using alkali-metal atoms. These have energies given by the Rydberg series E n * = 
— l/2n*, where n* = n — 5(n,l), and 5(n,l) is a quantum defect. The empirical 
parameters 5(n, I) give the energy-level shifts away from hydrogenic values. For large 
n, they approach asymptotic values 5(1) independent of n. 



In ref. [12], it is shown that the RSS construction can be generalized to include the 



effects of quantum defects. The analysis uses a model called supersymmetry-based 
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quantum-defect theory (SQDT) to describe alkali-metal atoms [30|. This analyti- 
cal theory for alkali-metal atoms has exact asymptotic quantum-defect energies as 
eigenvalues. The SQDT eigenfunctions form a complete and orthogonal set. 

The SQDT eigenstates are R n *i*(r)Yi m (9, 4>), where n* = n—5(l), I* = l—5(l)+I(l), 
and 1(1) is an integer that depends on I. The radial eigenstates R n *i*(r) have the same 
functional form as the hydrogenic functions R n i(r), but n is replaced by n* and I is 
replaced by I*. The angular wave functions are the usual eigenstates Yi m (6, 0) of the 
angular momentum. 

Our KSS construction in three dimensions can be generalized to include the effects 



of quantum defects pif . Since the angular part of the solution in SQDT separates, 
the SSS wave functions remain unchanged and are given by (ft) in Eq. (|I3D . The 
RSS wave functions ip(r) are obtained by writing the classical SQDT hamiltonian in 
terms of new radial oscillator variables incorporating effects of quantum defects and 
then finding minimum-uncertainty solutions for the corresponding quantum problem. 
The resulting wave function ip(r) has a related functional form to that in Eq. (p2|) 



and is discussed in ref. HT2 



A KSS ^f(r,9,4>) for alkali- metal atoms is again formed as a product of an RSS 
ip(r) and an SSS The functional form of the solution is similar to that in 

Eq. ( |2"3"D but with suitable replacements for the quantum numbers. To allow for the 
shifted energy eigenvalues of the alkali-metal atoms, the initialization procedure for 
the parameters must be modified. We choose (p r ) = and (L 3 ) = (3, and specify AL 3 
as before. Denote by E n * the energy expectation of the wave packet in the excited 
alkali-metal atom, and let r* ut be the outer apsidal point for a superposition of states 
with quantum-defect eigenenergies. Then, we impose (H) = E n * and (r) = r* ut , which 
differ from the hydrogenic case. To calculate (H) explicitly, we can take advantage 
of the completeness of the SQDT eigenfunctions and expand the initial KSS as a 
superposition of SQDT eigenstates, 

*(r,e,(j>,0)=^c nlrn Il n . l .(r)Y lm (e,(f>) , (32) 

n,l,m 

where the expansion coefficients, which depend on the KSS parameters, can be de- 
termined by inversion. The expectation value for the hamiltonian is then specified 
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by 

(H) = ^2 \c n i m \ 2 E n * = E n * . (33) 

n,l,m 

It is known that the long-term revival times for an alkali-metal wave packet depend 
on the quantum defects and that the effects of the quantum defects are different from 



the effects of a laser detuning po| . In addition, the appearance of deviations from 
the hydrogenic potential arising from the presence of core electrons in an alkali-metal 
atom means that the Runge-Lenz operator A does not commute with the hamiltonian. 
The classical orbit therefore precesses at a rate determined by the quantum defect 

0- 



V. SUMMARY 

In this paper, we have obtained minimum-uncertainty wave-packet solutions for 
the Schrodinger equation with a Coulomb potential in three dimensions. The solu- 
tions are the KSS and are given as a product of RSS with SSS. The RSS, previously 
derived, minimize the uncertainty relation for radial variables expressing the radial 
Coulomb problem in oscillator form. The SSS, constructed here, minimize the uncer- 
tainty relation for angular-coordinate and angular-momentum operators. The KSS 
provide analytical solutions to the Coulomb problem that move along classical kep- 
lerian orbits. They exhibit both classical and quantum-mechanical features. 

The KSS can be used as an analytical tool for studying the quantum-classical 
correspondence in the Coulomb problem. They may also be used to describe Ryd- 
berg wave packets created by excitation of a Rydberg atom with a short laser pulse 
in the presence of external fields. Such wave packets are expected to move in three 
dimensions along elliptical orbits that are strongly peaked around a plane. To match 
a KSS to a Rydberg wave packet, we choose the outer apsidal point as the initial 
location of the wave packet. The five KSS parameters are determined from the ex- 
pectation values of the radial position, the radial momentum, the energy, the angular 
momentum L3 transverse to the orbital plane, and the spread in L3. 

We obtained the time evolution of a KSS and provided an explicit example. The 
wave packet moves along an elliptical orbit with the classical keplerian orbital period. 
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The width of the wave packet oscillates during the motion, as is characteristic of 
a squeezed state. The KSS maintain their shape for several orbital cycles before 
collapsing and undergoing quantum recurrences. 

Finally, we provided an extension of the construction to the case where quantum 
defects are present. With this analysis, the KSS can be used for the description of wave 
packets in alkali-metal atoms, which are the ones of choice in current experiments. 
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APPENDIX 

In this appendix, we discuss the functions Aj(S) defined in Eq. (fUD and examine 
some of their properties. 

The integral in Eq. (|TJ|) can be evaluated. We find 

^W=4v^E^(f)r(A; + i)^ +fc (25) , (34) 

where i n (z) is a modified spherical Bessel function. This exact expression permits the 
numerical computation of Aj(S) to arbitrary precision using standard procedures. 

A useful identity for the A? (5) can be derived by taking advantage of some prop- 
erties of the modified spherical Bessel functions. We obtain 

Ap\8) = A$ +2 (S) + 3 -±±A§ +1 (6) . (35) 

For j > 0, we find A^(0) = 0. For j = 0, however, the value of the function with zero 
argument is 

p _ 4.(2/3)!! 
Ao{0) - (2/3+1)!! • (36) 
These relations are used to simplify some expressions in the main body of the text. 
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FIGURE CAPTIONS 



Fig. 1: A slice through the initial KSS wave packet in the x-z plane, which is trans- 
verse to the plane of the classical orbit. The quantity r 2 \^(r, 9, <p)\ 2 (in arbitrary 
units) is plotted as a function of x and z at t — 0. 

Fig. 2: Slices through the KSS wave packet in the x-y plane, at different times during 
the classical orbital cycle. The quantity r 2 |\l/(r, 9, <f), t)\ 2 (in arbitrary units) is 
shown as a function of x and y at the times (a) t — 0, (b) t — |T c i, (c) t = |T c i, 
(d) t = |T cl , (e) t = T ch 

Fig. 3: The KSS wave packet at time t = |T cl shown in Fig. 2c, viewed from a point 
on the positive x axis looking toward the nucleus. 
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